A direct calculation of critical exponents of two-dimensional anisotropic Ising model 
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Using an exact solution of the one-dimensional (ID) quantum transverse-field Ising model (TFIM), 
we calculate the critical exponents of the two-dimensional (2D) anisotropic classical Ising model 
(IM). We verify that the exponents are the same as those of isotropic classical IM. Our approach 
provides an alternative means of obtaining and verifying these well-known results. 
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It is well known that the two-dimensional (2D) classi- 
cal Ising model (IM)* has a thermodynamic continuous 
phase transition (CPT) from (anti-)ferromagnetic phase 
to paramagnetic phase at a critical temperature. Later it 
was also known 2 *— that a one-dimensional (ID) quantum 
transverse-field Ising model (TFIM) can be mapped into 
the 2D anisotropic classical IM under certain limit. Thus, 
a connection between the ground state of ID quantum 
TFIM and thermodynamical equilibrium state of a 2D 
classical IM was established. So far, our knowledge about 
the thermodynamical properties of 2D IM, especially 
those about critical phenomena, were mainly obtained by 
some sophisticate techniques such as high-temperature 
and low-temperature expansion—, renormalization group 
calcualtion 5 , Onsager's exact matrix-calculation and its 
variations 6 ' 7 . Normally, these techniques are computa- 
tional involving, and it is often not easy to be understood 
by the beginners and people outside of the field. In this 
short paper, we would like to introduce a direct way to 
obtain the quantities like critical exponents of 2D IM. 
The method is based on the exact solution of the ground 
state of ID TFIM 8 , which, in turn, can be obtained by 
using a simple Jordan- Winger transformation—. 

The Hamiltonian of a 2D classical IM is 1 
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where a XiV is the classical Ising spin on site (x,y) with 
two possible values ±1, and J x and J v are spin-spin cou- 
plings along x— axis(row) and y— axis (column) directions, 
respectively. It is well known that this model has a ther- 
modynamic CPT at a critical temperature T = T^£. 
Later it was shown 2 '- that with the relations and limits 



J; x 

krTf 



k B T 



1 - 2 Jy 

= — p k R T 



+oo 





(2) 



the transfer matrix between two neighboring columns 
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where 1 is identity operator and 
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is the Hamiltonian operator of a ID quantum TFIM 
where and a n are standard Pauli operators on site 
n. Thus the 2D anisotropic classical IM is mapped into 
a ID quantum TFIM with the following correspondences 
between the two models* : (1) the equilibrium state of 
2D IM vs. the ground state of ID TFIM; (2) the free en- 
ergy of 2D IM vs. the ground-state energy of ID TFIM; 
(3) the ensemble averages of physical quantities of 2D IM 
vs. the ground-state expectation values of time-ordered 
operators of ID TFIM; (4) the temperature T of 2D IM 
vs. the transverse field h of ID TFIM. 

With the so-called hard-core boson transformation 
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the Jordan- Wigner transformation— 
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the Hamiltonian of ID TFIM is rewritten into 



H = Nh 



dkHk 
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where N is site number, c n ,c n ,Ck,c\., are fermion opera- 
tors. Denote |0) as fermions' vacuum state, the ground 
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state and ground-state energy of H are obtained as£ 
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the excitation gap of H is 

A B = 2|J||t-l| 
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Thus the critical point of the ground-state CPT of the 
ID TFIM is at t = t c = 1 where Ae vanishes. 

Now that the ground state of ID quantum TFIM is 
obtained, critical behaviors of physical quantities of the 
ground state can be calculated directly. Let us first con- 
sider the 'heat capacitance'. According to the correspon- 
dences between 2D IM and ID TFIM, the 'heat capaci- 
tance' in ID TFIM is 
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Put in the ground-state energy obtained in Eq.® and 
JTJ3J, we have 
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where 



g(t) = (t-t c )/Vt 
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Near the critical point, i.e., t « t c = 1, the denominator 
in the integral can be expanded as 
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After some calculations one can obtain 
/ 1 32 
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Therefore, the critical behavior of the heat capacitance 
is the same as that of 2D isotropic classical IM 6 . 



The magnetization along the spin-spin interaction di- 
rection in ID TFIM is defined as 
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However, there are two facts which make it unable to ob- 
tain the correct value of M x directly. One is that the ID 
TFIM in Eq.|@J has no magnetic field component along 
the spin-spin interaction direction so that M x defined in 
Eci. (|18|) is always zero. The other is that in fermion 
representation is a non-local operator so that it is very 
difficult to calculate its expectation value even when the 
ground state of ID TFIM with non-zero magnetic field 
component along the spin-spin interaction direction is ob- 
tained. Fortunately, the magnetization can be extracted 
from the ground-state spin-spin correlation function de- 
fined as& 
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This correlation contains two parts of contributions. The 
first part is a long-range correlation of constant value due 
to the long-range magnetic order of the ground state, and 
the amplitude of this part is the square of the magneti- 
zation. The second part is the correlation of quantum 
fluctuations, and the amplitude of this part goes to zero 
in the limit r — > oo. Thus the magnetization equals to 
the square root of the amplitude of the correlation func- 
tion in Eq. (|19fl in the limit r — > oo. The spin-spin cor- 
relation function of ID TFIM was evaluated extensively 
in Barouch and McCoy's paper—, where the correlation 
function was expressed in the form of a Toplitz deter- 
minant and its asymptotic behaviors at large yet finite 
distance r were calculated with the use of a theorem of 
Szego. One exact result is 8 
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where sign(J) is the sign of J. Here we provide a direct 
way to obtain this result. Let us consider a finite chain 
of N spins with free boundary condition and calculate 
the correlation function of the first spin and the N— th 
spin, i.e., the ground-state expectation value of <rf <7^ in 
the limit TV — > oo. With Eq.Q, Eq.® and Eq.0, it is 
easy to show 
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Thus the rest work is to transform (cj + &i)(c N — cjv) 
into k— space representation by Eq.Q and calculate its 
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ground-state expectation value. From the form of |$g > 
in Eq. © , it is obvious that only the following four kinds 
of terms have non-zero expectation values 
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After doing some integrals, one finally obtains that in the 
thermodynamic limit N — > oo 

< QeffidfrlQa >= [sign{J)] N {l - i 2 ) 1/4 (24) 

for t < t c = 1 and is zero for t > t c , which is the same as 
Ea. (|20|l . From this result, one can see that 

M x = (1 - t 2 ) 1 /* « (t c - t) 1 /* (25) 

for t c — t — ► +0 and is zero for t > t c , which is exactly 
the same as the critical behavior of the magnetization of 
2D isotropic IMS. Thus we have found that the critical 
behaviors of the heat capacitance and the magnetization 
of ID TFIM and 2D classical isotropic IM are the same. 
According to the scaling laws of critical phenomena—, we 
can conclude that all the six critical exponents of the two 
models are equal. Thus our approach serves as an alter- 
native way to obtain and verify the well-known results of 
2D isotropic IM. In fact, another exact result obtained by 
Barouch and McCoy^, the value at the transition point 
t = t c which decays as a power-low function of r, 

G xx (t = t c ,r) = [sign(J)] r r~ 1 l A e 1/A 2 1 / 12 A-3 oc r" 1 / 4 
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where A = 1.282427130 is Glaisher's constant, was shown 
to be the same as that of 2D isotropic classical IM. The 



exponential decay of the quantum-fluctuation-induccd 
spin-spin correlation around the transition point was also 
shown to be equivalent to that in 2D isotropic IM&. It 
should be noted that it is unable to obtain the critical 
behaviors of zero-field susceptibility x an d the magneti- 
zation at the critical point under vanishing field because 
the ID TFIM in Eq.QJ has no magnetic field component 
along the spin-spin interaction direction. 

We would like to make some discussions about the 
above results. From the mapping between 2D classical 
IM and ID TFIM, one can see that the approximation in 
Eq.© is valid only when the couplings along one direc- 
tion, say, J y , goes to positive infinity while the couplings 
along the other direction, say, J x , goes to zero. This 
means that a ID TFIM is mapped into a 2D classical 
IM with extremely anisotropic couplings. Furthermore, a 
ID TFIM with anti-ferromagnetic coupling, i.e., J < 0, 
should correspond to a 2D classical IM with infinitely 
large ferromagnetic coupling along one direction and van- 
ishing anti-ferromagnetic coupling along the other direc- 
tion. Thus the equivalence of the critical behaviors of 
ID TFIM and 2D classical isotropic IM shows that the 
critical behaviors of 2D classical IM is very robust. 

In summary, with the use of the exact solution of the 
ID TFIM, we calculate the critical exponents of the 2D 
anisotropic classical IM and verify that the exponents 
are the same as those of 2D isotropic classical IM. This 
shows that the universality class of critical behaviors in 
2D classical IM is quite robust, since the 2D classical IM 
corresponding to a ID TFIM is extremely anisotropic. 
Our approach provides an alternative means to obtain 
and verify those well-known results. 
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